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Abstract 

Fundamental to many transportation network studies, traffic flow models can be 
used to describe traffic dynamics determined by drivers' car-following, lane-changing, 
merging, and diverging behaviors. In this study, we develop a deterministic queueing 
model of network traffic flow, in which traffic on each link is considered as a queue. 
In the link queue model, the demand and supply of a queue are defined based on the 
link's fundamental diagram, and its in- and out-fluxes are computed from junction flux 
functions corresponding to macroscopic merging and diverging rules. We demonstrate 
that the model is well defined and discuss its relationships with existing link-based 
and kinematic wave models. With numerical simulations of traffic dynamics in a 
diverge-merge network, we discuss the the consistency between the link queue model 
and the kinematic wave model. As an addition to the multi-scale modeling framework 
of network traffic flow, the model strikes a balance between mathematical tractability 
and physical realism and can be used for analyzing traffic dynamics, developing traffic 
operation strategies, and studying drivers' route choice and other behaviors in large-scale 
road networks. 

Keywords: Link queue model, link-based models, kinematic wave models, fundamental 
diagram, link demand and supply, junction flux functions, macroscopic merging and diverging 
rules 

1 Introduction 

Fundamental to many transportation network studies, traffic flow models can be used to 
describe traffic dynamics determined by drivers' car-following, lane-changing, merging, and 
diverging behaviors, subject to constraints in network infrastructure and control measures. 
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For a network traffic flow model, its inputs include initial traffic conditions, traffic control 
signals, and traffic demands determined by drivers' choice behaviors in routes, destinations, 
departure times, modes, and trips, and its outputs include vehicles' trajectories and travel 
times as well as congestion patterns. 

A traffic flow system, which is highly complex due to heterogeneous and stochastic 
characteristics of and interactions among drivers, road networks, and control measures, can 
be modeled at different spatio-temporal scales: at the vehicle level, microscopic models 



have been proposed to describe movements of individual vehicles (Gazis et al. , 1961 Gipps 



1986; Nagel and Schreckenberg , 1992 Hidas, 2005); at the cell level, kinematic wave models 



(Lighthill and Whitham, 1955; Richards, 1956), hydrodynamic models (Payne, 1971 Whitham 



1974), and gas kinetic models (Prigogine and Herman, 1971) have been proposed to describe 



the evolution of densities, speeds, and flow-rates inside road segments; at the link level, 
models based on variational formulations (Newell, 1993 Daganzo 2006), exit flow functions 
(Merchant and Nemhauser, 1978 Friesz et al. 1993 Astarita, 1996), and delay functions 



Friesz et al. , 1993) have been proposed to describe the evolution of traffic volumes on 



individual links; and at the regional level, the two- fluid model (Herman and Prigogine 



1979) and macroscopic fundamental diagrams (Daganzo and Geroliminis, 2008 Geroliminis 



and Daganzo! 2008) have been proposed for static traffic characteristics, and continuous 



models have been proposed to describe the dynamical evolution of the traffic density on a 
two-dimensional plane (Beckmann, 1952 Ho and Wong 2006). These models have different 
levels of detail and form a multi-scale modeling framework of network traffic flow: different 
models can describe different traffic phenomena at different spatio-temporal scales, and 
models at a coarser scale are usually consistent with those at a finer scale on average. 

Traditional vehicle- and cell-based traffic flow models have been widely applied in studies 



Daganzo 


1996; 


Lo 


1999) 



detailed to be mathematically tractable for many transportation problems in large-scale road 
networks, such as dynamic traffic assignment problems (Peeta and Ziliaskopoulos , 2001). In 
contrast, existing link-based network loading models are more amenable to mathematical 
analyses but fail to capture critical interactions among different traffic streams when queues 
spill back at oversaturated intersections (Daganzo, 1995a). In this study, we attempt to 



fill the gap between kinematic wave models and network loading models by proposing a 
link-based deterministic queueing model. The model is consistent with both kinematic wave 
models and existing link-based models but strikes a balance between mathematical tractability 
and physical realism for network traffic flow. Here we consider traffic on each link as a queue, 
and the state of a queue is either its density (the number of vehicles per unit length) on a 
normal link or the number of vehicles on an origin link. Based on the fundamental diagram 
of the link, we define the demand (maximum sending flow) and supply (maximum receving 
flow) of a queue. Then the out-fluxes of upstream queues and in-fluxes of downstream queues 
at a network junction are determined by macroscopic merging and diverging rules, which 
were first introduced in network kinematic wave models. Hereafter we refer to this model as 
a link queue model. 
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It is well known that a traffic system can be approximated by a deterministic network 
queueing system, in which traffic dynamics are dictated by link characteristics and interactions 
among traffic streams at merging, diverging, and other bottlenecks (Newell, 1982). In the 
transportation literature, point queue models have been used to model traffic dynamics on a 



road link (Vickrey, 1969 Drissi-Kai'touni and Hameda-Benchekroun, 1992; Kuwahara and 



Akamatsu 



1997). Point queue models are similar to fluid queue models for dam processes 



proposed in 1950s (Kulkarni, 1997). For stochastic queueing models of network traffic flow, 



refer to (Osorio et al. 2011) and references therein. Different from existing queueing models, 
the link queue model is deterministic, link-based, and highly related to kinematic wave 
models by incorporating the fundamental diagrams and macroscopic merging and diverging 
rules of the latter. The model captures important capacity constraints imposed by links 
and junctions but ignores the detailed dynamics on individual links. Therefore, the model is 
suitable for analyzing traffic dynamics, developing traffic operation strategies, and studying 
drivers' route choice and other behaviors in large-scale road networks. 

The rest of the paper is organized as follows. In Section 2, we review link-based models 
and kinematic wave models of network traffic flow. In Section 3, we present the link queue 
model and discuss its analytical properties and a numerical discrete form. In Section 4, we 
examine the relationships between the model and existing models. In Section 5, we apply 
it to simulate traffic dynamics in a simple road network and compare the model with the 
kinematic wave model. In Section 6, we discuss future research topics. 



2 Review of network traffic flow models 

For a general road network, e.g., a grid network shown in Figure [TJ the sets of unidirectional 
links and junctions are denoted by A and J, respectively. If link a G A is upstream to a 
junction j G J, we denote a — > j; if link a is downstream to a junction j, we denote j — > a. 
The set of upstream links of junction j is denoted by A^j = {a & A \ a — > j}, and the set of 
downstream links of junction j is denoted by Aj_> = {a G A | j — > a}. If a ^ A/_> for any 
j G J, link a is an origin link; if a ^ A^j for any j G J, link a is a destination link. We denote 
the sets of origin and destination links by O and R, respectively. We have that O fl R = 0, 
O G A, and R C A. We denote the set of normal links by A', where A' = A \ {O U R). For 
a G A', its length is denoted by L a . Here origin and destination links are dummy links with 
no physical lengths. 

In a traffic system, vehicles can be categorized into commodities based on their attributes 
such as destinations, paths, classes, etc. For examples, in dynamic traffic assignment problems 
without considering heterogeneous drivers and vehicles, vehicles on the same path belong 
to the same commodity, and a path flow is equivalent to a commodity flow; in emergency 
evacuation situations, vehicles who comply with evacuation orders may not have pre-defined 
routes, and the commodity of such vehicles can use many paths; if multi-class vehicles or 
multiple lane-groups are considered, then vehicles of two commodities can share the same 
path. A network with only one commodity is called a single-commodity network, in which 
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Figure 1: A grid network 



vehicles have identical characteristics and no predefined destinations or paths; otherwise, it is 
called a multi-commodity network. The set of commodities in the whole network is denoted 
by Q. If commodity w G fi uses link a G A, we denote u ~ a. The set of commodities using 
link a is denoted by fi a ; i.e., Q a = {u G Q \ u ~ a}. Then a unidirectional traffic network 
can be characterized by A = (A, O, R, J, {(A^j, Aj-+) : j G J}, {fl a '■ a £ ^4})- 

2.1 A link-based modeling framework based on traffic conserva- 
tion 

For a traffic network A, we denote the average density on a normal link a G A' at time t 
by k a {t) and denote the in- flux and out- flux of link a by f a (t) and g a (t), respectively. For 
commodity u> on link a, we denote its average density by k atUJ (t), in- flux by f a ,uj(t), and 
out-flux by g a ,w{t)- Thus we have for a G A': k a {t) = X^ e n a k a ,w(t), f a (t) = XLea* W> 
and g a (t) = J2iven a d^wif)- For an origin link o G O, we denote the queue length at time 
t by -Ko(^) and denote the arrival rate (in-flux) and departure rate (out-flux) by f D (t) and 
g (t), respectively. For commodity uj on link o, we denote its queue length by K QtUJ (t), 
in-flux by f ,uj(t), and out-flux by g a ,uj(t). Thus we have for o G O: K {t) = XLeo D Ko&{t)i 
fo(t) = Eo;en / ,^(*)) and 9o{t) = E^e n go,uj{t)- In this sense, we have a multi-commodity 
queue on each link: for each commodity, the number of vehicles on the link is the queue 
length, which equals the density times the link length on a normal link, the in-flux is the 
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arrival rate, and the out-flux is the departure rate. These link queues would interact with 
each other at junctions and result in a queueing network. 

From traffic conservation on normal and origin links, we have the following dynamical 
system for a traffic network A: (a G A' and o G O) 

dkJf) ' (fa(t)-9a(t)), (la) 



dt L a 
dk a Jt) 1 



dt L 
dK (t) 
dt 

dK , u {t) 



(fa ~ 9a A^) > (lb) 

fo(t)- 9o (t), (lc) 



(1/ foAt)-9oAt)- (Id) 

Further, from traffic conservation at a junction j, we have 

E = E m*)> ( 2a ) 

E 9a At) = E M*)- ( 2b ) 

Here we assume that there is no queue on a destination link r G R, and its in-flux equals the 
out-flux all the time. Then Q and ([2]) constitute a general link-based model of the queueing 
network A. It is a finite-dimensional dynamical system, whose dimension equals the total 
number of commodities on all origin and normal links. The link-based model, (JIJ and ([2]), is 
only based on traffic conservation on all links and junctions, and all traffic flow models on 
a network A should satisfy these conditions, whether they are vehicle-, cell-, or link-based 
models. 

However, the in-fluxes and out-fluxes, except the origin arrival rates, in ([T]) and ^ are 
under-determined. Therefore, ([TJ and ^ are incomplete, and various traffic flow models 
have been proposed to complement Q and §2§. The additional components are derived, 
explicitly or implicitly, from vehicles' car-following, lane-changing, merging, and divering 
behaviors on links and at junctions. In a complete traffic flow model, boundary fluxes and, 
therefore, queue lengths can be uniquely determined from initial conditions and boundary 
conditions in origin arrival rates and destination constraints. Thus traffic flow models have 
been used to predict traffic conditions on the whole network, which are essential for a better 
design, management, and control of a road network. 



2.2 A review of link-based network loading models 

In traffic assignment problems on a network A, all vehicles have pre-defined paths, and 
vehicles on different paths belong to different commodities. We denote the cumulative arrival 
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and departure flows of commodity u in the whole network by F u (t) and G w (t), respectively. 
Since commodity u has only one origin o (u 6 !]„ with o G O) and one destination r (u> G fl r 
with r G R), we have 

F w (*) 
G w (t) 

We denote the initial number of vehicles of commodity u in the whole network by N u (0), which 
equals N u (0) = K OjU (0) + J2 aeA' a)6n„ ^o,w(0)-^o- If we assume that vehicles of commodity ui 
follow the first-in-first-out (FIFO) principle^] then the travel time of a commodity u vehicle 
entering the network at t, IFj(t), can be found from the following optimization problem: 

= minvr, s.t. G u (t + tt) = F u (t) + JV w (0). (3) 

Note that + tt) = + N u (0) could have multiple solutions in n when G w (t) is not 
strictly increasing. When G u {t) is strictly increasing, ^ is equivalen to 

n.(t) = Gz\F u (t)+N u (0))-t. (4) 

In this sense, any traffic flow model (TFM) satisfying the commodity or path FIFO rule 
constitutes a mapping from origin arrival rates and initial conditions into commodity traffic 
times: 



/ f ,u>(r)dT, to G o G O 
Jo 

/ 9r,Lu(^)dT, u>^Vl r ,rER 
Jo 



II(.) = TFM(f(-),N(0)), (5) 

where II(-) is the set of all commodity travel times, f(-) the set of all origin arrival rates, and 
N(0) the set of initial queue lengths on all links. Link-based traffic flow models, Q and 
([2j, satisfying the path FIFO principle are usually called network loading models for traffic 
assignment problems, in which selfish drivers choose their paths to minimize their own travel 



times (Wardrop, 1952). In the following we review existing network loading models, which 



differ from each other in the ways of complementing ([T| and (]2j. 

Static link performance functions. The link performance functions in the static 
traffic assignment problem can be interpretated as a network loading model ([5]) as 
follows. We assume that the origin arrival rates and traffic queue lengths on all links 
are constant during a peak period, and there are no origin queues. Then from ([!]) we 
have that k a>U) is constant for any a and to and that f a>U} = g a ^\ from (J5]) we have that 



1 Note that the commodity FIFO principle is a necessary condition of the link FIFO principle, but not 
sufficient. That is, the commodity FIFO principle is a weaker requirement. 
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fa,u = fo,u if w ~ a and u ~ o. Thus we have 

fa = 2 /o^, a E A' 



(6a) 
(6b) 



In the traditional static traffic assignment problem, the link travel time function is 
defined as 



n„ 



(6c) 



and the travel times of all commodities are the same on a link. In addition, the link 



travel time Il a in (6c) can also be defined as a function of total flow during the peak 



period. Furthermore, the commodity travel time is 



(6d) 



Therefore, in the static traffic flow model ([6]), Q and (|2j) are not needed, and we can 
directly calculate link time-independent in- and out-fluxes and the commodity travel 
times from the assignments f 0jUI at origins. From the link travel time functions r a (/ a ), 
we have the following observations on the static traffic flow model (|6j): (i) traffic is 
static, (ii) link queues do not interact with each other, since the travel time on a 
link is irrelevant to flows on adjacent links; (iii) there are no capacity constraints on 
in- and out-fluxes; and (iii) vehicles follow the FIFO principle on all links, since all 
commodities have the same travel times on the same link. Various extensions have been 
proposed in the literature to address many limitations of the static traffic flow model, 
and the corresponding traffic assignment problems become asymmetric, multi-class, and 



capacitated (Boyce et al. 2005). Among the extensions are dynamic network loading 



models in the following. 



Delay functions. In (Friesz et al. , 1993), a delay function model was proposed to 



complement (jl]) and (|2j). On link a e A', one can obtain a relation similar to ([3]): 

F a (t) + L a k a (0) = G a {t + U a {t)), 

where F a (t) = Jq f a ( T )dT and G a (t) = g a (r)dr are cumulative link in- and out-flows, 
respectively. The corresponding differential equation is 



fa{t) = 9a(t + U a (t))(l + tl a {t)). 

Further, a dynamic link performance function is assumed: 

n (t) = r a (k a (t)). 



(7a) 



(7b) 
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The two equations above, coupled with (|TJ and Q, yield a link-based network loading 
model. In this model, there are no origin queues, and, at origin o G O, we have 
9o,ui(t) = fo,uj(t)- More discussions on theoretical and numerical properties of the model 



and its limitations can be found in (Daganzo, 1995a; Nie and Zhang 2005b Carey and 



Ge, 2007) and references therein. Generally, network loading models based on delay 
functions have the following features: (i) traffic dynamics and travel times on a link are 
purely determined by the in-flux and the initial conditions; (ii) all commodities have 
the same travel times on the same link; and (iii) in- and out-fluxes are not bounded by 
link capacities. 



Exit flow functions. In (Merchant and Nemhauser, 1978), an exit flow function was 
introduced as follows: 

9a(t) = Ga(k a (t)). (8) 

In the network loading model based on ([T]), (|2]), and Q, origin queues are not considered, 
and one can determine g a {t) and k a (t) from f a {t) and calculate path travel times from 

a 



Carey 


1986 


Friesz et al. , 


1989 



and references therein. This model share the same features as the one based on ([7j). In 



(Carey, 2004), some extensions of exit flow functions were proposed to incorporate queue 
spillbacks, but they are limited without considering merging and diverging behaviors 
at network junctions. 



• 


Point queue models. In ( 


Vickrey 


1969 


Drissi-Kaitouni and Hameda-Benchekroun 




1992; 


Kuwahara and Akamatsu, 


1997 


), vehicles are assumed to always travel at the 



free-flow speed on a link but wait at the downstream end before leaving the link. 
Therefore, at any time t, vehicles on link a can be divided into two groups: one group 
travels on the whole link at the free- flow speed V a , and the other group waits at the 
downstream end. In this sense, a link queue can be separated into a point queue and a 
traveling queue. We denote the density of traveling vehicles by k' a (t) and the number 
of waiting vehicles by X a (t). Thus we have 

X a (t) = L a -{k a {t)-k' a {t)). (9a) 

For the traveling vehicles, we have 

dk' a (t) 



dt 



^(fa(t)-fa(t-^)), 

i-ln V n. 



since all vehicles travel at the constant speed, V a . If we set the initial k' a (0) 
and A o (0) = 1^} Then the out-flux of link a is determined by 

min{/ a (t - C a (t)}, when X a (t) = 
C a (t), when X a (t) > 



9a(t) 



(9b) 

fca(0), 

(9c) 



2 Usually it is assumed that fc a (0) = 0, but impacts of initial conditions on the model's well-definedness 
are to be examined. 
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where C a (t) is the maximum discharging flow-rate of the link. Generally C a (t) is set as 
the time- independent capacity of the link (Smith, 1993). Then we have 



dXgjt) 
dt 



fa(t 



Va 



9a{t). 



(9d) 



Thus Q and Q yield the point queue model. Compared with network loading models 
based on delay functions and exit flow functions, the point queue model is more realistic, 
since it can capture capacity constraints on out-fluxes. But the out- flux and the point 
queue length of a link queue are determined by upstream in- fluxes, not downstream 
constraints caused by queue spillbacks. In (Nie and Zhang, 2005a), a spatial queue 



model was proposed to partially remedy queue spillback issues, but these models are 
still limited without fully considering merging and diverging conflicts among queues at 
network junctions. 



Even though congestion effects are modeled in existing network loading models, inter- 
actions among merging and diverging traffic streams are not fully captured, and vehicles' 
movements inside a link are not explicitly modeled. In addition, it is challenging to calibrate 
link performance functions, exit flow functions, or maximum exit flow-rate functions in these 
models. Due to insufficient physical realism in these models, they have only been applied 
in studying dynamic traffic assignment problems and not in analyzing traffic dynamics or 
developing traffic operations strategies. 



2.3 A review of network kinematic wave models 

In network kinematic wave models, which are extensions of the LWR model to network traffic 
flow, traffic dynamics inside a link can be described by the evolution of traffic densities at all 
locations. For link a G A' in a network A, we can introduce link coordinates x a G [X~,X£], 
where X~ and are upstream and downstream boundaries, respectively. Here the positive 
direction of x a is the same as traffic direction, X~ < X+, X~ > — oo, and X+ < oo. Then 
a road network with link coordinates can be written as A x n a eA'KTj -^a1> ^ n which any 
location can be uniquely determined by the link coordinate (a, x a ). For a G A', at a point 
(a,x a ) and time t, we denote the total density, speed, and flow-rate by p a (x a ,t), v a (x a ,t), 
and q a (x a ,t), respectively; we denote density, speed, and flow-rate of commodity uj G Vl a 
by PaA x a,t), v a>UJ (x a ,t), and q a ,u(xa,t), respectively. Here < p a (x a ,t) < p a ,j(x a ), where 
p a ,j(x a ) is the jam density at x a . 

Then a link-basedj^Jmulti-commodity kinematic wave model can be defined by the following 
rules with (x a ,t) suppressed from all variables unless necessary: 

(Rl) Additive densities and flow-rates: p a = XLen a and <?a = XLef^ 

3 We can also have a commodity-based model when different path flows belong to different commodities. 



9 



(R2) The constitutional law of commodity and total flows: q a ^ = p a ,uVa,u), and q a = p a v a . 
From Rl and R2, v a is the space-mean speed of all commodity speeds. 

(R3) Local and instantaneous fundamental diagrams: speed-density relations v a)U3 (x a ,t) = 
V a)U (x a , t,p a (x a ,t)), where p a is the vector of all commodity densities. 

(R4) Traffic conservation: %^ + ^ = 0, and ^ + §^ = 0. 



From R1-R4, we can obtain the following commodity-based LWR model on a normal link 
a G A': 



dpa,uj(^X a ,t^j dp au) (x a , t) Va,ui \^ai 1 ) ? ") ) 



dt 



dx a 



0. 



(10) 



which is a system of hyperbolic conservation laws defined on A x FJ A [X~ , X£]. Similarly 
we can obtain a hyperbolic conservation law on each link in total density, which is the 
traditional LWR model. But for link a there are only P a = |Q a | indepedent equations. 



The multi-commodity kinematic wave model, (10), encompasses many existing models: the 
homogeneous LWR model for single-commodity traffic on a link with a location- and time- 
indepedent fundamental diagram; the inhomogeneous LWR model for a single-commodity 



traffic on a link with location-dependent fundamental diagrams (Lebacque, 1996; Jin and 



Zhang, 2003a); the multi-lane-group model for heterogeneous traffic on a link (Daganzo 



Zhang 


2002 


); 


2002 


); 


2012a) 



; the general multi-class model for heterogeneous traffic on a link (Wong and Wong 



1997 



; and the multi-commodity network model for single-class traffic on a network (Jin 



However, discontinuous shock waves can develop in hyperbolic conservation laws even 
with continuous initial data, and there can exist multiple solutions (jSmoiler 1983). Therefore, 



in order for the network kinematic wave model (10) to be well-defined, we have to introduce 
two more rules: 



(R5) Weak solutions: (10) can admit non-differentiable shock wave solutions. 



(R6) Entropy conditions: some physical laws have to be used to pick out unique, physical 
weak solutions. 

Among the six rules, (Rl), (R2), (R4), and (R5) are the same for all systems of hyperbolic 
conservation laws in different application areas, but the fundamental diagrams in (R3) and 
entropy conditions in (R6) have to be specifically developed for network traffic flow. Generally 
the two rules can be derived and calibrated from static and dynamic characteristics of various 
traffic systems. 

Here we focus on the multi-commodity kinematic wave model of single-class, single-lane- 



group network traffic flow developed in (Jin, 2012a). That is, vehicles of different commodities 



have the same characteristics, and there is a single lane-group on each link. Therefore vehicles 
at the same location share the same speed with a speed-density relation of v au = v a = 
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V a (x a ,p a ). The corresponding flow-density relation is q a = Q a (x a ,p a ) = PaV a (x a ,p a ), and 
Qa, w = ^a,ujQai where commodity tu s proportion is ^ a ,uj = Pa,w 

jp a . Generally, Q a (x a ,p a ) is a 
unimodal function in p a and reaches its capacity, C a (x a ), when traffic density equals the 



critical density p a ,c{x a ) ( Greenshields 1935 Del Castillo and Benitez, 1995). Then (10) can 
be simplified as 



•^ai Pa) 

dt dx a 
dp a dp a V a ( Xaj Pa ) 



0. 



0. 



(11a) 
lib) 



dt dx a 

At (x a ,t), in addition to traffic density, speed, and flow-rate, two new variables, traffic 



demand and supply, can be defined as follows (Daganzo, 1995b): 



d a 



D a (x a ,p a ) = Q a (x a ,mm{p a , p atC (x a )}), 

Sa\ X ai Pa) = Qai^ay^^-iPa: Pa,c( x a)}), 



(12a) 
(12b) 



where the demand d a increases in p a , and the supply s a decreases in p a . For single-class and 
single-lane-group traffic flow, commodity demands are proportional to commodity densities; 
i-e., d a ,oj(x ai t) = d a (x ai t) pa p '^ a ^ . Then at a junction j, out-fluxes, gj(t), and in-fluxes, fj(t), 
can be computed from upstream commodity demands, dj(t), and downstream supplies, Sj(t), 
using the following flux function: 



MMt)) = ¥F(d s (t),sj(t)). 



(13) 



Note that (13) can also be used to determine fluxes through a location inside a link. These 



flux functions were first introduced in the Cell Transmission Model (CTM) to describe 



priority-based merging and FIFO diverging rules (Daganzo, 1995b), and many other flux 



functions have been developed in the literature for other macroscopic merging and diverging 
behaviors. In ( |Jin et al.[ [2009] [jln] |2010b|a] |2012a[ ), it was shown that ([131 constitutes a valid 
entropy condition (R6) in the sense that the Riemann problem for (11) is uniquely solved 



within the framework of (Holden and Risebro, 1995). In a sense CTM can be considered a 



first-order discrete form of the network kinematic wave model, (10) and (13). In addition, 



traffic dynamics at origins can be treated as follows. At origin o, if there is no queue, its 
commodity demands are the arrival rates; otherwise, the commodity demands are set to be 
very large. That is, we define an origin link's demand by 



d {t) 



oo, K (t) > 
f (t), K o (t) = 



(14) 



and (13) can still be used to determine fluxes through the junction downstream to the origin. 



Therefore the network kinematic wave model completes (jTJ and (J2j) with detailed kinematic 
waves on all links. 

In summary, the kinematic wave model of network traffic flow has the following features: 
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• A link's storage capacity is determined by the jam density, the in- flux of a downstream 
link does not exceed its local supply, and the out-flux of an upstream link does not 
exceed its local demand. 

• On a link, vehicles may follow the FIFO principle or not. 

• Speed-density relations on a link and merging and diverging rules are explicitly defined 
in demand and supply functions and junction flux functions. 

• For a road network with given origin demands, destination supplies, and initial condi- 
tions, the network kinematic wave model can be understood as an infinite dimensional 
dynamical system, in which state variables are traffic densities at all locations. 

• Although the network kinematic wave model can describe detailed wave propagations 
in a road network, it does not capture drivers' delayed responses, hysteresis in speed- 
density relations, or other properties of microscopic car-following models. That is, 
they are continuous approximations of car-following models by retaining important 
characteristics on links and at network junctions. 

The network kinematic wave model has been used to study traffic dynamics, operations, or 



traffic assignment problems in road networks (Daganzo 1996 Jin, 2009). It is an important 



component in the multi-scale modeling framework of traffic flow. However, being an infinite- 
dimensional dynamical system, it is both computationally and analytically demanding. 



3 A link queue model 

In this section, we present a new link-based model to complement Q and ([2]). We consider 
traffic on a link as a single queue and call this model as a link queue model. For each 
link queue, the state variable is the link density, k a (t) (a e A'), or the link volume, K a (t) 
(o G O). In this model, in-fluxes, out-fluxes, and travel times can be computed from link 
densities. Different from existing link-based models, this model inherits two major features 
from network kinematic wave models: first, the local fundamental diagram is used to define 
the demand and supply of a link queue; second, flux functions are used to determine in- and 
out-fluxes from link demands and supplies at all junctions. But different from the network 
kinematic wave model, the new model only describes the average traffic dynamics on a link 
and does not capture shock or rarefaction waves. 

3.1 The link queue model for single-class, single-lane-group traffic 

In this subsection, we consider single-class, single-lane-group network traffic and further 
assume that a normal link a (a e A') is homogeneous^] with a local fundamental diagram 

4 An inhomogneous road can be divided into a number of homogeneous ones. 
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q a = Qa(Pa) at all locations for p a G [0,fc QJ ], where k a j is the jam density on link a. In 
addition, Q a (Pa) is a unimodal function in p a , and the capacity C a is attained at critical 
density k afi \ i.e., C a = Q a {k ac ) > Q a (p a )- F° r a normal link queue a (a G A'), we extend the 
definitions of local demands and supplies in (12) and define its demand by 

Qa(k a (t)), k a (t) G [0,k a , c ] 

C a -, k a (t) G {k a>c ,k a j] 



daif) 

and its supply by 

Sa(t) 



<3 a (min{fc a (t),fc a)C }) 



(15a) 



Q a (max{£; a (t), fc a , c }) 



(15b) 



C a , A; a (t) G [0, fc a , c ] 

Qa{k a {t))i k a (t) G {k a ,ci k a j] 

For commodity w, its proportion is denoted by £ 0)W (t) = k a ^(t) / 'k a (t) , and its demand is 
proportional to ^ a ,w(^)- 

For an origin link o G O, if we omit the origin queue, then its demand, d a (t), and the 
commodity proportions, ^ 0)W (t), should be given as boundary conditions. Otherwise, if the 
arrival rates / (t) and f ,u(t) at the origin are given as boundary conditions, a point queue 
can develop at the origin, and we define its demand by 

K {t) > 



doit) = fo(t) + I Ko 



(*)>o 



oo, 







(15c) 



where the indicator function lK a (t)>o is infinity if K a (t) > and zero otherwise. For a 
destination link r G R, its supply, s r (t), is given as boundary conditions: if the destination 
link is not blocked, we can set s r (t) = oo. 

At a junction j, we apply (13) to calculate corresponding in- and out-fluxes from upstream 
demands and downstream supplies 

fe(t),fi(t)) = FF(d,(t),s,(t)), (16) 

where dj(t) is the set of upstream commodity demands, Sj(t) the set of downstream supplies, 
gj(t) the set of out-fluxes from all upstream links, and fj(t) the set of in- fluxes to all 
downstream links. 

Therefore, completing ([TJ by demand-supply functions in (15) and well-defined flux 
functions in (16), we obtain the following link queue model of network traffic flow (a £ A' 
and o G O): 

1 (U B (k(f)) - 3.(k(f))) , (17a) 



dfcq(t) 
dt 

dk atUJ (t) 

dt 
dK (t) 
dt 

dK^Q 
dt 



L a 

1 



(U a , w (k(i))-5V(k(t))), 



f (t)-Q (k(t)), 

foAt)-3oAHt)), 



(17b) 

(17c) 

(17d) 
(17e) 
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where k(i) is the set of all link densities or volumes, and f a {t) = U a (k(£), g a (t) = 3 a (k(i)), 
fa^tf) = ^(k(t)), g a>u {t) = 9 a ,a>(k(t)), = 3 (k(t)), and flr , w (*) = 9 0(W (k(*)) are 



computed from k(t) with (jl5j and (16). The link queue model, ( |17[ ), is a system of first-order, 
nonlinear ordinary differential equations, and the number of state variables equals the number 
of commodities on all normal and origin links, J2 a eA' \^a\ + Soeo When the initial 

states and boundary conditions are given, state variables at all times can be calculated from 



3.2 Some examples of junction flux functions 



A flux function (16) should be consistent with physically meangingful merging and diverging 



rules. A well-defined flux function in (16) should have the following properties: 



1. Traffic conservation at a junction, ([2]), is automatically satisfied. 

2. A link's out-flux is not greater than its demand. As a special case, if a link's demand is 
zero, its out-flux is zero. 

3. A link's in-flux is not greater than its supply. As a special case, if a link's supply is 
zero, its in-flux is zero. 

4. The flux function should be Godunov or invariant for the network kinematic wave 



model in the sense of (Jin, 2012b). That is, non-invariant flux functions, which can be 



used in network kinematic wave models, cannot be used in the link queue model. 

In this subsection, we present some valid flux functions for four types of typical junctions 
shown in Figure [2} These flux functions have been proved to be invariant or Godunov for 



Jin et al. 


2009 


Jin 


2010b|a 


2012alb 



(16) can be defined as follows: 



9i{t) 



f 2 (t)=mm{d 1 (t),s 2 (t)}, 



(18a) 

;isb) 



Here we assume that the boundary fluxes are maximal subject to constraints in the upstream 
queue's demand and the downstream queue's supply. When signals or other types of 
bottlenecks present at the junction, this model has to be modified. 

For a merge with two upstream links, 1 and 2, and a downstream link, 3, as shown in 
Figure [2](b), if the two upstream links are the same type of roads, e.g., freeways, (16) can be 
defined by the fair merging rule: 



f 3 {t) = min{di(t) + d 2 (t),S3(t)}, 
gi(t) = min{<ii(t), max{s 3 (t) — d 2 (t), 



Ci 



C X + C 2 



ss(t)}}, 



(19a) 
(19b) 
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(a) A linear junction 



(b) A merge 





(c) A diverge (d) A general junction 

Figure 2: Four types of junctions 



92(t) 



fs{t)-gi{t), 

2 



(19c) 
(19d) 



a=l 



Note that the non- invariant fair merge model in (Jin and Zhang 



cannot be used in (16). If the two upstream links have different merging priorities; e.g., when 



2003b), g x (t) = min{di(t), s 3 (t)^ 



di(t) 



(t)+d 2 (t)- 



an on-ramp has a higher merging priority than a freeway, (16) can be defined as follows: 

f 3 (t) = wm{di(t) + d2(t),s 3 (t)}, (20a) 
gx(t) = min{di(t),max{s 3 (t) - d 2 (t),as 3 (t)}}, (20b) 
92(f) = fs(t)- gi (t), (20c) 



(20d) 



(1=1 



where a is the merging priority of link 1. Obviously the fair merging model is a special case 
of the priority merging model. 

For a diverge with an upstream link, 0, and two downstream links, 1 and 2, as shown in 
Figure |2^c), if all vehicles have pre-defined route choices and follow the FIFO diverging rule, 
(16) can be defined as follows: 



g (t) = minH(t),-^-,-^-}, 



(21a) 
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hit) = £o-i(*)(7o(*), (21b) 
hit) = Mt)go{t), (21c) 
f a ,Ut) = 9oit)£oAt), a = 1,2, u G tt a ,u G tt . (21d) 

where f -n(*) = Eueflonn! and £ ^ 2 (0 = Eu,en nn 2 6,^(0 are the proportions of 

vehicles on link traveling to links 1 and 2, respectively. Here £ ->i(0 > 0, £o->2(£) > 0, 
and £ ->i(£) + £o^2(£) = 1- Again, the non-invariant diverge model in (Lebacque, 1996), 
flit) = min{£o-»i(i)rfo(£), cannot be used in the link queue model. In emergency 

evacuation situations, vehicles have no pre-defined route choices, (16) can be defined as 
follows: 

g (t) = mm{d (t),s 1 (t)+s 2 (t)}, (22a) 
Mt) = mm{ Sl it),max{do(t)-s 2 it),f3d it)}}, (22b) 
/ 2 (t) = goit)-fi(t), (22c) 

where /3 is the evacuation priority of link 1. 

For a general junction j with m upstream links and n downstream links, as shown in 



Figure |2Jd), if all vehicles follow the FIFO diverging and fair merging rules, (16) can be 
defined as follows (A->j = {1, • • • , m} and = {m + 1, • • • , m + n}): 

1. From total and commodity densities on all links, k a (t) (a G A->j), k au ,it) (a G A_^, 
uj G f2 a ), fcb(t) (6 G A(_j), and k^^it) (6 G A^-, u; G fife), from (15) we can calculate 
all upstream demands, <i a (t) (a G A_>j), downstream supplies s&(t) (6 G A/-*-), and the 
turning proportion £ ->b(i) 

= M*)> ( 23& ) 

where EbeA- = 1 for a G ^4^j. Note that origin demands and destination supplies 

could be given as boundary conditions. 

2. The out-flux of upstream link a G A_»j is 

g a (t) = min{4(t),^(t)C a }, (23b) 
where the critical demand level 9 jit) uniquely solves the following min-max problem 

^•(t) = mm{ max {^^j, mm max — — }. (23c) 

Here A\ it) a non-empty subset of A_» 3 - . 

3. The commodity-flux is (a G A_>j, 6 G Aj_>., a; G fi a fl 

/mW = g a ,Ut) =9ait)UUt)- (23d) 
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4. The in- flux of downstream link b G is 



h(t) 



9a(t)£a->b(t). 



aeA_ 



(23e) 



Since (18), (19), and (21) are its special cases, (23) is a unified junction flux function. 
We have the following observations on the unified junction model (p3j): (i) From (23e), 



in- and out- fluxes satisfy the conservation equations (j2|; (ii) The model satisfies the first 
in- first-out (FIFO) diverging rule (Daganzo, 1995b), since the out-fluxes of an upstream 
link are proportional to the turning proportions; (iii) The model satisfies the fair merging 
rule (Jin, 2010b), since, when all upstream links are congested, d a {t) = C a for a G A 
9a{t) = Oj{t)C a < C a , and the total out-flux of link a is proportional to its capacity 



3.3 Properties 



In this subsection we focus on the link queue model defined by (17), (15), and (23), where 



vehicles have predefined routes and follow the fair merging and FIFO diverging rules. Then 
we obtain the following finite-dimensional link queue model: 



dk(t) 
dt 



F(k(t),u(t);0), 



(24) 



where u(t) denote boundary conditions in origin demands or arrival rates and destination 
supplies, and include link lengths, fundamental diagrams, speed limits, metering rates, 
numbers of lanes, and other network and driver characteristics. In an extreme case with 



triangular fundamental diagrams, if all links carry free flow, then (24) becomes a linear 



system and is therefore well-defined. Here we demonstrate that (24) is well-defined with (15) 



and (23) under general traffic conditions. 



Lemma 3.1 The critical demand level in (23c) is well-defined: its solution exists and is 



unique. In addition, 0j(t) G [0,1], and it is a continuous function of upstream demands, 
downstream supplies, and turning proportions. 



Proof. Since the number of non-empty subsets of A^j is finite, the min-max problem in (23c) 



has a unique solution. Thus Oj(t) has a unique solution for any combinations of upstream 
demands, downstream supplies, and turning proportions. Obviously it is a continuous function 

in these variables. In addition, 8j(t) > 0, since m&x^^cA ^ ' 

s b {t) 



i_kj(t) C<»£a->(>(*) 



> 0. Since d a (t) < C a from (15), we have 9j(t) G [0, 1]. 



> 



Lemma 3.2 In (23), < g a (t) < d a (t), and < fb(t) < Sb(t). That is, in- and out-fluxes 



are bounded by the corresponding demands and supplies. 
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Proof. From (23 
As shown in 



Jin 



we hav e < g a (t) < d a (t). 

2012b], 6j{t) = max a6 ^. ^ if and only if s b (t) > £ ae . 



for all b; in this case, g a (t) = d a (t), and f h (t) = £ oeA d a (t)£ _>. 6 (t) < s b (t) for all b. 



Otherwise, there exists A* such that 



0j(t) = 
9 3 {t) < 
9 ^t) > 



mm 

dg(t) 

c a ' 

dg(t) 



a G A* 



a G A^j \ A* 



Then from (23b), we have g a (t) = 9j(t)C a for a G A* and g a (t) = d a (t) for a G A^j \ A*. 



Thus for all b G A~ 



hit) 



c ^ b (t) + d ^)^ b (t) < Sb (t). 



Therefore < f b (t) < s b (t). ■ 

Lemma 3.3 g a is a continuous function of d a . In addition, it is piece-wise differentiable in 
d a : 



dg a 
dd a 



1, d a < 9jC a 
0, d a > 9jC a 



Therefore, g a is Lipschitz continuous in d a . In addition, all out- and in-fluxes in (23) are 
Lipschitz continuous in all upstream demands and downstream supplies. 

Proof. Since 9j is continuous in d a , g a is also continuous in d a . From properties of 9j discussed 
in Section 4.3 of (Jin, 2012a), we can see that: (i) 9j increases in d a when d a < 9jC a ; (ii) when 

d a > 9jC a , 9j is constant at 9* = min^^. — T^~^o~f — z — > where a G A\. Therefore, 



when d a < 9*C a , from (23b) we have g a = d a , and = 1; when d a > 9*C a , from (23b) we 
have g a = 9*C a , and ^ = 0. Similarly, we can use properties of 9j to prove that all out- and 
in-fluxes in (23) are Lipschitz continuous in all upstream demands and downstream supplies. 



Theorem 3.4 The link queue model of network traffic flow, (17) together with (15) and 
(23), is well-defined. That is, under any given initial and boundary conditions, solutions to 
the system of ordinary differential equations (24) exist and are unique. 
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Proof. For general flow-density relations in fundamental diagrams, both traffic demand and 



supply defined in (15) are Lipschitz continuous in traffic density. Further from Lemma (3.3), 



we can see that the right-hand side of (24) is Lipschitz continuous in all state variables 



as well as boundary conditions. Then from the Picard-Lindelof theorem, solutions to the 



system of ordinary differential equations ( 24 ) exist and are unique under any given initial 



and boundary conditions (Coddington and Levinson, 1972). That is, the link queue model is 
well-defined. ■ 



In summary, the link queue model, (24), has the following properties: 



1. If a link is empty, from (15) its demand is zero, and its out-flux is zero. Thus from (17) 



the link's density is always non-negative. 



2. If a link is totally jammed, from (15) its supply is zero, and its in-flux equals zero. Thus 



from (17), the link's density cannot increase after it reaches the jam density. Therefore, 



a normal link's density is always bounded by its jam density. 

In addition, a link's out- and in-fluxes are bounded by its demand and supply, re- 
spectively. Since demand and supply are not greater than the capacity, the in- and 
out-fluxes are also bounded by road capacities. 

At a junction, all related link queues interact with each other. Especially when a 
downstream link is congested, even if it is not totally jammed, the upstream links will 
be impacted due to the limited supply provided by the downstream link. Therefore, 
queue spillbacks are automatically captured. 

The information propagation speed on a link may not be finite, since traffic on a link is 
always stationary instantaneously in a sense. Thus the model cannot capture shock or 
rarefaction waves inside a link. 



For link a, the dynamic link flow-rate q a can be defined as q a = Q a {k a ). But it is not 
used in the model, and may not be the same as the in-flux and out-flux, even when 

traffic is stationary; i.e., when ^jj^ = 0. Thus link flow-rates are less important than 
in- and out-fluxes in the model. Similarly, the link travel speed q a /k a is not explicitly 
considered. 



3.4 A numerical solution method 



The link queue model, (24), cannot be analytically solved under general initial and boundary 
conditions, but many numerical methods are available for finding its approximate solutions 
(Zwillinger, 1998). Here we present an explicit Euler method for the model. We discretize 
the simulation time duration [0, T] into M time steps with a time-step size of At. At time 
step i, the total and commodity densities on link a G A' are denoted by k l a and k l aLU ', K l Q and 
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K l ou denote the numbers of vehicles, i.e., queue lengths. In addition, the boundary fluxes 
during [iAt, (i + l)At] are denoted by /*, /* g\, and g* . 



On a normal link a, its demand, d l a , and supply, s*, can be computed from k l a with (15). 
For an origin link o, its demand can be computed as follows: 

< = (25) 



Then traffic states at time step i + 1 can be updated with the discrete version of (24) 



= K + P(fa-9% (26a) 

C = k U + P- (fi»-9iJ, (26b) 
K +1 = K + (r o -gl)At, (26c) 
= iC,. + (/;.-< J At, (26d) 

(26e) 



where the boundary fluxes are computed by (23) with densities at time step i; i.e., at a 
junction j, the boundary fluxes for all upstream links a G A^j and downstream links b G 
are given by 

*U = E I (27a) 



ft* 



6 - = mm max — , mm max — j—^. k 27b 

£ = mm{di,ep a }, (27c) 

£ = E (27d) 

= 9a,u = 9a(a,u) (27e) 

where A\ a non-empty subset of A_> 3 -. Here the arrival rates /* and /* w and destination 
supplies are given as boundary conditions. Note that, for an explicit Euler method, the 



time step-size At should be small enough for the discrete model, (26), to converge to the 
continuous version. In addition, the smaller At, the closer are the numerical solutions to 
theoretical ones. 



20 



4 Relationships with existing models 



4.1 Relationships with existing link-based models 



The link queue model, (17), can be considered as an extension to the exit flow function model 
(J8J), since out-fluxes are computed from link densities. However, in the link queue model, 
in-fluxes are also calculated from link densities, and both in- and out-fluxes are determined 
by densities of all links around a junction. 

In the link queue model, a point queue model is used for an origin link. But different 



from the traditional point queue model, we define the demand in (15c), and the out-flux is 



determined by the flux function (13) at the junction downstream to the origin queue. That is, 
the origin out-flux is determined by downstream links' supplies when they're congested. But 
in the traditional point queue model, interactions between a point queue and its downstream 
queues are not fully captured. 

In summary, compared with existing link-based network loading models, the link queue 
model has the following properties: 

• The densities are bounded by jam-densities. That is, k a (t) G [0, k a j] on a normal link a. 

• Speed-density and flow-density relations are directly incorporated into the demand and 
supply functions. 

• At a junction, merging and diverging rules are included, and interactions among different 
links are explicitly captured. 



The FIFO principle is automatically satisfied in such junction flux functions as (23). 



Link travel times can be calculated from in- and out-fluxes but are not included in the 
model. 



4.2 Relationships with the network kinematic wave model 



The link queue model, (17), can be considered an approximation of the network kinematic 



wave model (11), since (i) fundamental diagrams of the kinematic wave model are used to 
calculate link demands and supplies and, therefore, in- and out-fluxes, and (ii) invariant flux 
functions of the kinematic wave model are used to calculate in- and out-fluxes through a 
junction. However, different from the kinematic wave model, the link queue model cannot 
describe the formation, propagation, and dissipation of shock and rarefaction waves on links, 
since traffic is instantaneous stationary on each link. 

Note that, in (Yperman et al. , 2006), a link transmission model was proposed based on 



the variational version of the LWR model by (Newell. 1993). In their model, the fundamental 



diagrams and merging and diverging rules are also consistent with those in the kinematic 
wave models, but shock and rarefaction waves are still captured, only implicitly. Thus 
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the link demand and supply functions in that study are different from the demand and 
supply functions of a link queue in this study, and their model can be considered as another 
formulation of network kinematic wave model and is still infinite-dimensional. 



The discrete equations in (26) are equivalent to the corresponding multi-commodity CTM 



when each link is discretized into only one cell. This also suggests the consistency between 
the link queue model and the kinematic wave modelj^] From this relationship, we conclude 
that t he time-step size, At, in (26) should satisfy the following CFL condition for all normal 
links (ICourant et al.L |l928b, Vj^< 1; i.e., 



At 



aeA> V a 



(2* 



where V a is the free-flow speed on link a. That is, the maximum time-step size should not be 
greater than the smallest link traversal time in a network. 



Following (Jin 



2003 Chapter 6), one can develop a simulation platform of network traffic 



flow based on ( 26 ) with appropriate network representations, data structures, and program 
flow-charts. Simulations can be carried out in such programming languages as C/C++ and 
Matlab. In addition, the continuous link queue model can be implemented with Matlab 
Simulink directly. For a road network with n origins and n destinations, the computational 
cost is in the order of n 3 / At, if there is one and only one path between an O-D pair [^J the 
computational cost for the corresponding CTM is in the order of n 3 / At 2 . In contrast, if the 
commodity densities are not tracked when turning proportions are known at all junctions or 



vehicles have no pre-defined routes, then the computational cost of (26) is in the order of 
/At[j 



n 



5 Numerical examples 



In this section, we apply the link queue model (24) to study traffic dynamics in a diverge- 
merge network with two intermediate links, referred to as the DM2 network, shown in Figure 
|3j In the network, there are two commodities: vehicles of commodity 1 use link 1; and those 
of commodity 2 use link 2. We denote the dummy link at the origin by 4 and that at the 
destination by 5. In this example, we do not consider the origin queue. Initially, the network 
is empty: k a = for a = 0, • • • , 3. Link lengths are L a — 1, 1, 2, 1 miles for a = 0, ■ ■ • , 3, 
respectively; the number of lanes are n a =3, 1, 2, 2, respectively; and all links share the 



following triangular fundamental diagram (Munjal et al. 



, 1971 


Haberman 


1977 


Newell 



5 In a sense, CTMs can be considered cell queue models. 

6 For example, for a grid network with n east-west and n north-south streets, the number of links on a 
path is proportional to n, and there are n 2 paths. Then the number of commodities densities is in the order 
of n 3 , and the computational cost is in the order of n 3 / At. 

7 Also consider a grid network with n east-west and n north-south streets, the number of links in the 
network is 2n(n + 1), and the computational cost is in the order of n 2 /At. 
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1 



3 




Figure 3: A diverge-merge network with one O-D pair and two intermediate links 



1993): 



lk k 
Qa(k a ) = mm{v f k a , -{n a - -^)} = min{65fc a , 2925(n a - t^t)}, 

T Kj 1 U 

where the free-flow speed Vf = 65 mph, the jam density kj = 180 vpmpl, and the time-gap 
r = 1/2925 h. Therefore, the critical density k atC = 36n a , and 

d a = min{65£; a , 2340n a }, 

s a = min{2925(n a -^),2340n a }. 



We solve the link queue model with the numerical method in Section 3^4 . The simulation 
time duration is T = 1.05 hrs, and At = 1.75 x 10 -4 hrs, for which the CFL condition (28) 



is satisfied. To compare the results with those of the kinematic wave model (11), we also 



solve the commodity-based CTM with the same fundamental diagrams, demand and supply 
functions, and merge and diverge models. For CTM, the cell size Ax = 0.0125 miles, and the 
corresponding CFL number = 0.91 < 1. 



5.1 Constant loading patterns 

In this subsection we compare the link queue model and the kinematic wave model with 
constant loading patterns. Here the boundary conditions are constant: the origin demand is 
constant <i 4 (t) = C = 7020 vph; the destination supply is also constant s 5 (t) = C 3 = 4680 
vph; and the proportion of commodity 1 at the origin is £ 4) i = £, where £ is constant but can 



take three different values: 0.3, 0.45, and 0.7. In (Jin, 2009), it was shown that the DM2 



network reaches damped periodic oscillatory, persistent periodic oscillatory, and stationary 
solutions respectively in these three cases. Since traffic dynamics are dictated by those 
on links 1 and 2 in the network, in the following we only present link densities^ in- and 
out-fluxes on these two links. 

In Figure [4| we demonstrate the results from the link queue model and the kinematic 
wave model when £ = 0.7 in the DM2 network. These figures confirm that, when £ = 0.7, 



3 In CTM, the link density equals the average value of all cell densities. 
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(a) (b) 




t (hour) t (hour) 

Figure 4: Comparison between the link queue and kinematic wave models when £ = 0.7: In 
all figures, the solid lines are results for the link queue model, and the dashed lines for the 
kinematic wave model; In figures (c) and (d), blue lines are for in-fluxes, and red lines for 
out-fluxes 
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traffic reaches stationary states on links 1 and 2 eventually, since traffic densities reach 
constant, and the in- and out-fluxes become equal on both links. We can observe the following 
similarities between the two models: 

• The two models have the same stationary states. 

• On average the two models share the same dynamical patterns in densities, in- and 
out-fluxes. 

• It takes a longer time for traffic to converge to stationary states on link 2 than on link 
1, since the former is longer. 

We can also observe the following significant differences between the two models: 

• Results from the link queue model converge in an exponential fashion, but those from 
the kinematic wave model converge in a finite time. 

• In the link queue model, traffic densities and fluxes on both links become positive 
immediately after traffic is loaded at t > 0; but in the kinematic wave model, it takes 
some time for traffic densities and fluxes on both links to become positive, since it takes 
time for vehicles to travel from the origin to the diverging junction. 

In Figure |5j we demonstrate the results from the two models when £ = 0.3. The dashed 
curves in figures (c) and (d) confirm that damped periodic oscillations occur on both links. 
But results from the link queue still converge to stationary states exponentially. In Figure 
[6j we demonstrate the results from the two models when £ = 0.45. The dashed curves in 
all figures confirm that persistent periodic oscillations occur on both links. The period is 
about 0.2 hours or 12 minutes, [j In both cases, results from the link queue still exponentially 
converge to stationary states, but the results from the link queue model are still consistent 
with those from the kinematic wave model on average. 



5.2 Varying loading patterns 

In this subsection, we compare the link queue model and the kinematic wave model with 
a varying loading pattern: the origin demand is periodic d±(t) = |Co(sin(47rt/T) + 1) vph, 
where T — 1.05 hrs; the destination supply is still constant s^(t) — C3 — 4680 vph; and the 
proportion of commodity 1 at the origin is £4^ = 0.45, which leads to persistent periodic 
oscillations with constant demands in the preceding subsection. Still, we only compare link 
densities, in- and out-fluxes on the two intermediate links. 

In Figure [7j we demonstrate the results from the two models. We can see that traffic 
dynamics on the network are dominated by the varying demand pattern. The dashed curve 



In (Jin 2009), it was shown that the period is determined by the lengths of links 1 and 2 as well as the 



corresponding fundamental diagrams. 
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Figure 5: Comparison between the link queue and kinematic wave models when £ = 0.3: In 
all figures, the solid lines are results for the link queue model, and the dashed lines for the 
kinematic wave model; In figures (c) and (d), blue lines are for in-fluxes, and red lines for 
out-fluxes 
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Figure 6: Comparison between the link queue and kinematic wave models when £ = 0.45: In 
all figures, the solid lines are results for the link queue model, and the dashed lines for the 
kinematic wave model; In figures (c) and (d), blue lines are for in-fluxes, and red lines for 
out-fluxes 
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Figure 7: Comparison between the link queue and kinematic wave models with varying 
demand patterns when £ = 0.45: In all figures, the solid lines are results for the link queue 
model, and the dashed lines for the kinematic wave model; In figures (c) and (d), blue lines 
are for in-fluxes, and red lines for out-fluxes 
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in all figures show that persistent periodic oscillations occur on both links when the traffic 
demand is higher than a certain level, and the period is still about 12 minutes. Clearly, even 
with varying demand patterns, the link queue model is still consistent with the kinematic 
wave model in the simulation results. 

From simulation results in this and the preceding subsections, we can see that the proposed 
link queue model approximates the kinematic wave model over a relatively long period of 
time, which can be in the order of 10 minutes. In particular, the two models are consistent 
in stationary traffic. Different from the kinematic wave model, in which traffic densities and 
fluxes can abruptly change over time, the link queue model yields much smoother quantities. 
In addition, as shown in Figure [6] and Figure [7j the kinematic wave model can be unstable 
in the DM2 network, but the link queue model is always stable. All these similarities and 
differences between the link queue model and the kinematic wave model of network traffic 
flow resemble those between car-following models and the LWR model of link traffic flow 10 
In this sense, the link queue model can be considered a continuum version of the kinematic 
wave model, since discontinuities caused by shock and rarefaction waves are smoothed in the 
former. 

6 Conclusions 

In this paper, we presented a link queue model of network traffic flow, in which the evolution 
of congestion levels on a road link is described by changes in the link density. With link 
demands and supplies, it can capture basic characteristics of link traffic flow, including 
capacity, free-flow speed, jam density, and so on. In addition, with appropriate junction flux 
functions, it can describe the initiation, propagation, and dissipation of traffic queues in a 
road network caused by merging, diverging, and other network bottlenecks. The link queue 
model is a finite-dimensional system of ordinary differential equations and highlights the 
most important two characteristics of network traffic flow: static fundamental diagrams and 
dynamic junction models. Thus, even though the model omits some details of link traffic 
dynamics, including shock and rarefaction waves, it retains the basic queueing properties 
of a link, including its holding capacity (jam density) and processing capacity (maximum 
flow-rate), and uses the merging and diverging rules to relocate the right of way at network 
bottlenecks. Both analytical and numerical results suggest that the model is a continuous 
approximation of the kinematic wave model in a large-scale road network during a time 
period in the order of 10 minutes. In summary, the link queue model strikes a balance 
between existing link-based network loading models and network kinematic wave models 
and is appropriate for many studies on congested large-scale networks, in which interactions 
among traffic dynamics on different links cannot be omitted. 

Even though the link queue model is more amenable to mathematical analyses, yields 

10 For example, on a ring road, some car-following models can be unstable, while the LWR model is always 
stable. However, they share the same steady states. 
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smoother travel time functions, and is computationally more efficient, it does not undermine 
the importance of kinematic wave models. In contrast, it highlights the importance of 



the latter, since the definitions of link demands and supplies in (15) and the junction flux 



functions in (23) are directly borrowed from existing kinematic wave models. We expect that 
studies on both types of models would benefit each other. Therefore, the link queue model is a 
useful addition to the multi-scale modeling framework of network traffic flow. In applications, 
we may first apply the link queue model to obtain analytical insights of network congestion 
patterns under different demand levels, control strategies, route choice behaviors, or other 
conditions and then apply the kinematic wave model as well as microscopic models to further 
study the propagation of traffic queues and other details before drawing any conclusions or 
making any policy recommendations. 

In the future we will be interested in developing link queue models of other traffic flow 
systems, which are consistent with kinematic wave models: 

• If commodity flows are not explicitly tracked, but the turning proportions at 
all junctions can be detected through loop detectors or other devices, we can obtain a 
link queue model of implicit multi-commodity traffic. In this case, only one equation, 



(la), is needed for the evolution of total traffic on a link; rt23p,c,e) can still be used to 



calculate in- and out-fluxes at a junction; and fb,u(t) and g a ,ui(t) are not available in 



(231). This model is suitable for traffic operations when route choice behaviors are not 



explicitly accounted for. 

If a network is closed without any origin or destination links, the link queue model can 
still be applied. In this case, turning proportions at all junctions can be exogeneous 
or endogenous, and the model becomes an autonomous system without boundary 
conditions in origin demands or destination supplies. 

The model can be extended for multi-class, multi-lane-group traffic systems with lane- 
changing traffic, HOV lanes, traffic signals, capacity drops, ramp metering, etc. The 



major challenge is to define traffic demands and supplies in (15) and extend the junction 



flux functions in (23) for such scenarios. 



With the link queue model, we will also be interested in studying the following problems 
pertaining to network traffic flow: (i) stationary states, or equilibria, of the link queue model 
in open or closed networks ( Jin |2012c ); (ii) the stability of the link queue model; (iii) hybrid 



link queue, kinematic wave, and car-following models; (iv) analyses and simulations of traffic 
dynamics in a large-scale road network with data input. In addition, the link queue model, 



(24), can be viewed as a control system, in which u are control variables. From the viewpoint 
of control systems, we can analyze the system's responses to control signals, and many 
transportation applications can be studied as control problems. 
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